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Abstract 

We present resistivity measurements of the upper critical field (ffc2) phase 
diagram as a function of temperature (T) for Ndi.85Ceo.i5Cu04+5/SrTi03 
single crystal films with different degree of disorder (6) in magnetic fields up 
to 90 kOe at temperatures down to 0.4 K. The data are well described by a 
two-band/two-gap model for a superconductor in the dirty limit. 
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1. INTRODUCTION 

The theory of the two-band superconductivity has been developed more 
than 50 years ago: Suhl, Matthias, and Walker [1] predicted the existence of 
multigap superconductivity, in which a distinction of the pairing interaction 
in different bands leads to different order parameters and to an enhancement 
of the critical temperature. 

A discovery of MgB2 has renewed interest in superconductors with mul- 
ticomponent order parameters. MgB2 seems to be the first superconductor 
for which a two-gap model offers a simple explanation of many anomalous 
experimental findings. The two-band model for superconductivity in MgB2 
was first suggested in [2,3]. On the basis of first-principles calculations of 
the electronic structure and the electron-phonon interaction, it was argued 
that superconductivity in this compound resides in two groups of bands: the 
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group of two strongly superconducting cr-bands and the group of two weakly 
superconducting vr-bands. 

A theory for the upper critical field Hc2 of two-gap superconductors in the 
dirty limit is developed by Gurevich [4] on the ground of a multiband BCS 
model [1]. The equations for Hc2 for two-band superconductors in the dirty 
limit are derived with the account of both intraband and interband scattering 
by nonmagnetic impurities. It is shown that the shape of the Hc2(T) curve 
essentially depends on the ratio of the intraband electron diffusivities Di and 
D2, and can be very different from the one-gap dirty-limit theory. Taking 
MgB2 as an example, the results of this work are used to describe the observed 
abnormal temperature dependences of Hc2 of this compound. 

The discovery of oxypnictides as a new class of superconductors [5] has 
regenerated interest in high temperature superconductivity. It is shown by 
Mazin et al. [6] that the pairing mechanism in iron arsenide pnictides is con- 
sistent with the multi-band theory: the so-called "extended s± -wave" model 
with a sign reversal of the order parameter between different sheets of Fermi 
surface. The relatively high upper critical fields (Hc2) with atypical temper- 
ature dependence in these materials immediately attracted much interest. 
In particular, the anomalous upward curvature observed in the temperature 
dependence of Hc2 in these materials made the dirty two-band model [4] used 
to describe MgB2 a natural and useful choice [7-11]. The model of Gurevich 
[4] is widely exploited for describing of a variety of the abnormal Hc2 tem- 
perature dependences, both for H || c and i/ ± c, in NdFeAsOo.yFo.s single 
crystals [7], in polycrystalline LaFeAsOo.sgFo.ii samples [8], in cobalt-doped 
SrFe2As2 epitaxial films [9], in Co-doped BaFe2As2 single crystals [10] and 
in Sri_a;Eu2;(Feo.89Coo.ii)2As2 single crystals (x = 0.20 and 0.46) [11]. 

An anomalous upward curvature of IIc2(T) dependence has been reported 
in the earliest investigations of electron-doped superconducting Nd2_2:Ce^Cu04_|_5 
crystals [12, 13]. Such a non-BCS behavior of Hc2(T) was observed later both 
in Nd2-a;Ce2:Cu04+5 single crystals [14] and thin films [15, 16]. 

Gantmacher et al. [14] have argued that the observed concave form of 
Hc2(T) in their bulk Ndi.82Ceo.i8CuO,i+5 crystals may be an indication of 
unconventional boson-type superconductivity. In the work of Golnik and 
Naito [15] it is shown that the normal state resistivity. Hall coefficient and 
magnetoresistivity of Nd2-a;Cea;Cu04+5 films can be described quantitatively 
within a simple two-carrier model if the existence of an electronlike and a 
holelike band is assumed. The physical origin of the two conduction bands 
at that time remained an open question. In the superconducting regime the 
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positive curvature of the resistive critical field was attributed by them to the 
quasi-2D character of the material and the resulting fiuctuation effects. 

There are much recent activities in investigation of possible electrons and 
holes coexistence in the normal state of cuprate superconductors [17]. Two 
kinds of carriers in electron-doped cuprates with different lanthanide cations 
seem to arise from the electronic structure near the Fermi surface (FS) of 
the Cu02 planes. Angle resolved photoemission spectroscopy (ARPES) has 
revealed a small electron-like FS pocket in the underdoped region, and a 
simultaneous presence of both electron- and hole-like pockets near optimal 
doping in Nd2-a;Cea;Cu04+5 systems [18,19]. The conclusions of the ARPES 
measurements and first-principle calculations of the electronic structure on 
the electron-doped high-Tc superconductors Lui 85Ceo.i5Cu04 (Ln = Nd, Sm 
and Eu) performed by Ikeda et al. [20] , are in accordance with these results. 
A spin density wave (SDW) model [21-23] was proposed which gives qual- 
itative explanation to ARPES observations. In this model, SDW ordering 
would induce FS reconstruction that results in an evolution from an electron 
pocket to the coexistence of electron-like and hole-like pockets and then into 
a single hole-like FS with increasing of doping. 

In our previous work [24] we have analyzed the Ce -doping dependence of 
the normal state Hall coeffient in optimally reduced Nd2-a;Cea;Cu04+5/SrTi03 
single crystal films and, on the grounds of this analysis, we have recruited 
a two-carrier model for describing of the mixed state Hall coeffcient. We 
have adopted a semi-phenomenological description of a mixed state Hall ef- 
fect by flux-flow model of Bardeen and Stephen modifled by coexistence 
of electrons and holes. In this work we present results on the upper criti- 
cal field as a function of temperature determined by a resistive method in 
Nd2_a;Ce2,Cu04+5/SrTi03 single crystal films with x= 0.15 and various de- 
gree of disorder (S). Our goal is to describe the observed abnormal Hc2 
(T) dependences on the ground of two-gap model of Gurevich [4]. Such an 
approach may occur actual in the light of much recent efforts on experimen- 
tal (ARPES) and theoretical (SDW-model) investigations of the electronic 
structure near Fermi surface in the Cu02 plane of electron-doped cuprate 
superconductors . 

2. SAMPLES AND EQUIPMENT 

The series of Nd2-a;Cea,.Cu04+5/SrTi03 epitaxial films {x = 0.15) with 
standard (001) orientation were synthesized by pulsed laser evaporation [16]. 
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The original target (the sintered ceramic tablet of Nd2-xCea;Cu04+5 of the 
given composition) was evaporated by a focused laser beam and the evap- 
orated target material was deposited on a heated single-crystal substrate. 
The substrate material was SrTiOs with (100) orientation and dimensions of 
5x10x1.5 mm. The substrate temperature was 800° C, the pressure during 
the deposition was 1.067 mbar, the residual gas was nitrous oxide (N2O). 
Then the films were subjected to heat treatment (annealing) under various 
conditions to obtain samples with various oxygen content. As a result, three 
types of samples with x=0.15 were obtained: as-grown samples, optimally 
reduced samples (optimally annealed in a vacuum at T = 780° C for t = 60 
min; p = 1.33-10"^ mbar) and non optimally reduced samples (annealed in 
a vacuum T = 780" C for t = 40 min; p = 1.33-10~^ mbar). Temperature 
dependence of resistivity were measured at the magnetic field up to 90 kOe 
using the Quantum Design PPMS in the temperature range T = (1.8 -j- 300) 
K and the Oxford Instruments superconducting magnet up to 120 kOe in the 
temperature range T — (0.4 4.2) K. 

3. EXPERIMENTAL RESULTS AND DISCUSSION 

The in-plane longitudinal resistivity pxx is measured as a function of tem- 
perature T for different applied magnetic fields H or as a function of H at 
fixed temperatures in single crystal films of electron-doped superconductor 
Nd2-a;Ce2;Cu04+5 in magnetic field perpendicular to the ab-plane up to 90 
kOe at T=(0.4 300)K. Fig. la shows the p(T) dependences for the op- 
timally reduced Ndi.85Ceo.i5Cu04 film. It is seen that an increase of the 
external magnetic field results in a shift of the resistive transition to lower 
temperatures without appreciable broadening. Thus the values of the upper 
critical field Hc2 can be deduced from p(T)) curves for fixed magnetic fields. 
For the definition of Hc2 we use a value that is the intersection of tangents to 
the resistive transition region and to the normal-state resistivity region (see 
Fig. lb). 

In Fig.2 we have plotted Hc2(T) dependences derived from measured p(T) 
and p(B) transitions for all our Ndi.85Ceo.i5Cu04_|_^ samples. In both opti- 
mally and nonoptimally reduced films we observe a positive curvature of 
Hc2(T) for temperatures from Tc down to about Tc/2 and then obvious ten- 
dency to saturation at low temperatures. 

On the insert of Fig.2 Hc2(T) dependences for optimally and nonoptimally 
reduced samples of Golnik and Naito [15] are presented for comparison. It is 
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Figure 1: The resistive transitions of the Ndi.85Ceo.i5Cu04/SrTi03 fihii with optimal 
anneahng: (a) as function of temperature at various magnetic fields; (b) as a function of 
magnetic field at various temperatures. 
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Figure 2: Temperature dependences of the upper critical field Hc2 defined as a value of 
the crossing of two tangents: first - to the resistive transition region; second - to the 
normal resistivity region. Inset: Hc2(T) dependences for optimally and nonoptimally 
reduced samples of Golnik and Naito [15] - open labels; our experimental data of Hc2(T) 
at T<4.2K - solid labels. 
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seen that their and our dependences have similar behavior at T > 4.2K and 
that our data rather well complement the data of [15] at T < 4.2 K. Mea- 
surements of the resistive critical field of Ndi.84Ceo.i6Cu04+5 single crystals 
have been reported for temperatures below 4.2 K with similar results (see 
[12,13]). 

As revealed by ARPES [18,19] the two types of carriers with opposite 
charges in the normal state of electron-doped cuprate comes from two parts 
of FS that is associated with a reconstruction of conduction band due to 
SDW ordering which splits the band into upper and lower branches [25, 26]. 
The inset of Fig. 3a shows the FSs for a typical optimally-doped sample 
[22]. The SDW correlation splits the continuum FS into two pieces of sheet: 
electron-like(l) and hole-like (2) pockets. 

In the SC state the reconstruction of Fermi surface leads naturally to a 
two-gap model. The situation with piecewise FS (with two types of FS parts) 
in superconducting state of electron doped cuprates is properly considered, 
for example, by Liu and Wu [22]. The model suggests two SC gaps with 
different amplitude: Ai for electron-like band and A2 for hole- like one. 
The two-gap model gives a unified explanation for the experimental data on 
ARPES [19] and Raman scattering [27] in the electron-doped cuprates with 
X < 0.17. 

On these grounds we explore theoretical considerations about tempera- 
ture dependence of the upper critical field in two-gap superconductor [4] for 
an interpretation of our experimental data. 

We regard Nd2-a;Cea;Cu04+5 single crystal with x= 0.15 as a supercon- 
ductor with two sheets (pockets) 1 and 2 of the Fermi surface on which the 
superconducting gaps take the values Ai and A2, respectively (indices 1 and 
2 correspond to electron and hole pockets) and the intraband diffusivities 
electron Di and hole D2. 

The equation for Hc2 for two-band/ two-gap superconductor in the dirty 
limit are derived in [4] with the account of both intraband and interband 
scattering by nonmagnetic impurities and paramagnetic effects. For a simpler 
case of negligible interband scattering and paramagnetic effects the equation 
for Hc2 has been presented as [4b]: 

ln{t) = -[U{h) + Uir]h) + A0H/2+ 

+ [iUih) - U{r^h) - \_/wf/A + \,2\2,/w'Y'\ (1) 

Here t — T/Tc, /i=Hc2Di/2$oT, $0 — 7r^/e - is a flux quantum, 77 = 
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D2/D1, 



[/(a:)=^(l/2 + x)-^(l/2), (2) 

where is a digamma function and X± — An ± A22, w— A11A22 - 

A12A21: Ao ={Xl + 4Ai2A2i)^/2^ 

Eq.(l)contains tlie matrix of the BCS couphng constants: the diagonal 
terms An and A22 describe intraband pairing (band 1 have the highest cou- 
phng constant An) and A12 and A21 describe interband couphng. 

For equal diffusivities, 1] = 1, Eq. (1) simplifies to the one-gap de-Gennes- 
Maki equation of Werthamer - Helfand - Hohenberg (WHH) theory (see [28] 
and references therein): 

ln(t) + U{h) = 0. (3) 

As was to be expected for Ai2=0 equation (1) breaks up on two indepen- 
dent WHH ones. For the first type of carriers with D=Di we have: 

/n(t(^)) + U{h) = 0, (4) 

where t^^^ =T/Ti^\ T^^^ ~ exp(-l/An). For the second type with D=D2 
we obtain from (1): 

ln{t^^^) + U{r]h)^\o/w, 

or 

ln{t^^^) + U{rih) = 0, (5) 

where t*^^) = T/Ti^-',Ti^-'~ exp(-l/A22) and expressions Ao/w=(Aii - A22)/(AiiA22) 
and T^^^/xi^^ = exp(-(Aii - A22)/(AiiA22)), rehable at A12 = 0, are used. 

On Figs 3, 4 we show the fit of the two-band theoretical equation (1) to 
Hc2(T) curves with fitting parameters listed in Table I for optimally (Fig. 
3a) and non optimally (Fig. 4a) reduced Ndi.85Ceo.i5Cu04+5 samples. We 
find that for both samples the curves of Hc2(T) can be rather well described 
by the two-band model for the diffusivity ratio rj = D2/Di< 1 and for nearly 
negligible interband BCS coupling constants: A12 <S An, A22- It follows from 
the previous transport measurements on two investigated samples [24] that 
carriers with the higher diffusivity Di are electrons and ones with the lower 
diffusivity D2 are holes. 
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Figure 3: Upper critical field versus temperature for optimally reduced Ndi.85Ceo.i5Cu04 
sample. The filled symbols correspond to the experimental data, (a) Solid line corresponds 
to Hc2(T) calculated from the two-gap Eq.(l) for the fitting parameters listed in Table 
1. (b) Red and blue curves are calculated according to the one-gap Eqs (4) and (5), 
respectively. The fitting parameters are listed in Table 2. 
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Figure 4: Upper critical field versus temperature for non optimally reduced 
Ndi.85Ceo.i5Cu04-(-5 sample. The filled symbols correspond to the experimental data, 
(a) Solid line corresponds to Hc2(T) calculated from the two-gap Eq.(l) for the fitting 
parameters listed in Table 1. (b) Red and blue curves are calculated according to the 
one-gap Eqs (4) and (5), respectively. The Aft ing parameters are listed in Table 2. 



Encouraging by the strong inequality between interband and intraband 
coupling constants we try to describe H^.2(T) dependences on the ground of 
independent equations (4) and (5) for A12 = 0. The fitting for optimally (Fig. 
3b) and non optimally (Fig. 4b) reduced samples shows that the observed 
dependences may be considered as a superposition of WHH-like curves for 
two types of carriers (electrons and holes) with different critical fields H^2^ 
or H^2^ and with different critical temperatures Tc^-* or The fitting 

parameters are presented in Table 2. As it is seen from Figs. 3b and 4b it is 
necessary to take into account the finite value of A12 only near a crossover of 
the two WHH curves at T = T cross- 

In accordance with analysis of Gurevich [4] in a case of different diffusiv- 
ities the limiting value of Hc2(0) is determined by the minimum diffusivity 
(i.e. by holes): 

Fe2(0) = $oTi^V27^2, (6) 

where 7 = 1.78. On the other hand the expressions for Hc2(T) near Tc 
is determined by the band with the highest coupling constant An (i.e. by 
electrons) : 

Hc2{T)=A^o{Tc-T)/7r^Di. (7) 
The experimentally observed ("real") values of Hc2 and Tc should be the 
highest ones: for critical fields these are Hc2^ ^ ^ T^ros-s and Hc2^ ^ 
> Tcross, for critical temperature Tc = (see Figs 3b, 4b). The lower 
("virtual" in terms of the article [4b]) values of Hc2 come to the hght, for 
example, in the magnetic field dependences of Hall coefficient, Rff(B), in a 
mixed state of superconductor (see [24]). 

On Fig.5 we present the data for Hc2(T) dependences in as grown Ndi 85Ceo.i5Cu04+5 
sample. In contrast to the other two samples we observe here a convex form 
of Hc2(T) curve at all temperature interval. The data are remarkably well 
described by the WHH-like equation in accordance with a result of transport 
measurements [24] that electrons and holes give nearly equal contributions to 
the conductivity in this sample, i.e. Di = D2. It is known [4] that for equal 
diffusivities of the two bands, i.e., rj = D2/D1 = 1, the parametric equation 
(1) reduces to the one- gap formula in the WHH theory. 

Prom an analysis of Rff(B) dependences in a mixed state [24] we have 
obtained that in as grown sample the critical field for electrons is higher than 
that for holes: li^^2 > '^^^ ■> both real (H^2'') ^^^d virtual (H^g'*) fields have 
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Figure 5: Upper critical field versus temperature for as grown Ndi.85Ceo.i5Cu04+5 sample. 
The filled symbols correspond to the experimental data. Red and blue curves are calculated 
according to the WHH formulas (4) and (5), respectively. The fitting parameters are listed 
in Table 2. 
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been determined (see Table 2). As for equal diffusivities one has T^^^/T^^^ 
= H^2VHc2'* ) ^6 can constructed not only real WHH curve for electrons but 
also virtual one for holes (see Fig. 5). 

It should be noted that the existence of two different spatial and mag- 
netic field scales, consistent with the experimental data on Abrikosov vortex 
structure in MgB2, has been demonstrated by Koshelev and Golubov [29]. 

They provided a quantitative model for the vortex structure in a two- 
band superconductor with weak interband impurity scattering and strong 
intraband scattering rates. It was shown that in a case of weak interband 
scattering the two typical sizes of the isolated vortex in different bands appear 
and, as a consequence, two magnetic field scales exist which can be accessible 
experimentally. 

4. CONCLUSIONS 

In conclusion, we report the data for the upper critical field as a func- 
tion of temperature for Ndi 85Ceo.i5Cu04+5 single crystal films with different 
degree of disorder (6) determined by the resistivity measurements. Two 
types of Hc2(T) dependences are observed: the abnormal ones with a pos- 
itive curvature of Hc2(T) for temperatures from Tc down to about Tc/2 in 
both optimally and non optimally reduced films and convetional WHH curve 
in as grown film. We find that the Hc2(T) dependences can be consistently 
explained by the two-band/two-gap model of a dirty superconductor for the 
diffusivity ratio r] = D2/D1 < 1 in the first two films and for equal diffusivi- 
ties, T) = D2/D1 = 1, in the third one. 

Thus, we have demonstrated that the two-band/two-gap model captures 
the essence of the behavior of Hc2(T) for investigated systems. This result 
is in accordance with experimental evidences and theoretical considerations 
for a coexistence of electrons and holes in a normal state of optimally doped 
Nd2-a;Ce3;Cu04+5 system. 

This work was done within RAS Program (project N 12-P-2-1018) with 
partial support of RFBR (grant N 12-02-00202). 
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Tabic 1: Parameters of the fits to the two-gap formula (1) for Ndi.85Ceo.i5Cu04-|-5/SrTi03 
epitaxial films. 



Samples 


T 


He2(0), 


( All Ai2\ 
VA21 A22/ 


Di 


D2 


ri 




K 


kOe 










Optimally 


20.2 


70.4 


/U.43 U.U44\ 
VO.044 0.35/ 


1.89 


0.45 


0.24 


reduced 














Non optomally 


17.3 


53.8 


/U.54 U.U38\ 
V0.038 0.44/ 


2.43 


0.57 


0.235 


reduced 
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Table 2: Parameters of the fits to the WHH-hke formulas (4) and (5)for 



Ndi.85Ceo.i5Cu04+5/SrTi03 epitaxial films. 



Samples 




Hi2^(0), 




Hg(0) 


Di 


D2 


V 




K 


kOe 


K 


kOe 








Optimally 


20.2 


30.0 


11.4 


70.4 


1.89 


0.45 


0.24 


reduced 
















Non optimally 


17.3 


20.0 


11.2 


53.8 


2.43 


0.57 


0.235 


reduced 
















As grown 


13.3 


12.2 


10.9 


10.0 


3.0 


3.0 


1 
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